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We consider the damped and driven dynamics of two interacting particles evolving in a symmetric 
and spatially periodic potential. The latter is exerted to a time-periodic modulation of its inclination. 
Our interest is twofold: Firstly we deal with the issue of chaotic motion in the higher-dimensional 
phase space. To this end a homoclinic Melnikov analysis is utilised assuring the presence of transverse 
homoclinic orbits and homoclinic bifurcations for weak coupling allowing also for the emergence of 
(3J0! hyperchaos. In contrast, we also prove that the time evolution of the two coupled particles attains 

^ , a completely synchronised (chaotic) state for strong enough coupling between them. The resulting 

'freezing of dimensionality' rules out the occurrence of hyperchaos. Secondly we address coherent 
collective particle transport provided by regular periodic motion. A subharmonic Melnikov analysis 
is utilised to investigate persistence of periodic orbits. For directed particle transport mediated by 
rotating periodic motion we present exact results regarding the collective character of the running 
solutions entailing the emergence of a current. We show that coordinated energy exchange between 
the particles takes place in such a manner that they are enabled to overcome — one particle followed 
by the other — consecutive barriers of the periodic potential resulting in collective directed motion. 
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I. LEAD PARAGRAPH 



The study of transport phenomena has attracted considerable interest over the years due to its 
relevance in many physical situations, which are often modelled on the basis of one-dimensional particle 
motion in a spatially periodic potential. If the particles, in addition to their motion in the periodic 
potential, interact, cooperative effects not found in situations of individual particle motion, may arise. 
The objective of the current work is to investigate the conditions under which it is possible to generate 
a directed flow along with collective motion in a system of two coupled particles. To be precise, we 
study the transport of two particles interacting via a harmonic bond force. The individual dynamics of 
uncoupled particle motion, evolving in a symmetric and spatiallay periodic potential whose inclination 
is time-periodically modulated, is chaotic, preventing the emergence of directed particle motion. For 
the dynamics of the coupled particles we elucidate the possible scenario in which the energy exchange 
between the particles proceeds in such a well-coordinated manner that the particles move separately 
| from one well into the next, one following the other, resulting in directed motion. Notably, we 
demonstrate, that the dynamics of the coupled particles exhibits coherent collective directed motion 
provided by transporting attractors. Furthermore, it is shown that for suitable coupling strengths 
between the particles, a completely synchronised chaotic state is attained confining the dynamics in 
phase space as a result of 'freezing of dimensionality'. 



II. INTRODUCTION 



Transport phenomena play a fundamental role in many physical systems. For a number of applications, including 



Josephson junctions charge density waves superionic conductors Q, rotation of dipoles in external fields 
phase-locked loops [5] and diffusion of dimers on surfaces |6]-|13f to quote a few, transport is based on the dynamics 
evolving in spatially periodic potential landscapes. When an additional external time-periodic modulation is applied 
to the periodic potential, interesting effects such as phase-locking, hysteresis [l4[ and stochastic resonance [l5] are 
found. Recent investigations have dealt with the Hamiltonian dynamics of individual particles evolving in a periodic 
potential whose inclination is time-periodically varied by a weak external monochromatic modulation field (l6l-tl8|. 
Remarkably, for the corresponding one-and-a-half degree of freedom Hamiltonian system it has been demonstrated 
that adiabatic modulations of the slope of the potential lead to the generation of transient transport dynamics related 
with enormous directed particle flow. An explanation for this behaviour has been given in terms of the underlying 
phase space structure of the externally driven one degree of freedom system promoting the motion in ballistic channels 
[19| . Recently interest in the collective transport dynamics of interacting particles in periodic potential landscapes 
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has grown (20j-(2^. 

In the present paper we consider the motion of two damped, time-periodically driven and coupled particles in 
a periodic potential. Our aim is to demonstrate that depending on the coupling strength between the two particle 
subsystems a rich dynamics results ranging from hyperchaos, stable periodic running states mediating coherent particle 
transport, to a completely synchronised chaotic state. Motion takes place in a five-dimensional phase space for which 
details of its intricate structures remain elusive not least due to the higher dimensionality. Whether in systems with 
a larger number of (microscopic) degrees of freedom such (macroscopic) behaviour as collective motion leading to 
a directed flow emanates from higher dimensional dynamics is not obvious as also more complex forms of unstable 
motion, namely hyperchaos may emerge psl ). We investigate the case when (already) the dynamics of a single 
(uncoupled) particle is characterised by deterministic chaos. In this context several questions arise: namely (i) if the 
chaotic motion persists in the coupled dynamics or, (ii) may the coupling drive the system's dynamics into a state of 
synchronisation, and (iii) may chaos even become suppressed (see |26|.|27|)? 

On the other hand, regarding coherent collective transport leading to the emergence of a directed flow it is paramount 
that a regime of regular periodic motion exists which is characterised by transporting attractors. Such a regime of 
regular periodic motion needs to emerge from the coupled dynamics (we recall that a single, uncoupled system 
exhibits chaos). These are the problems which we address. The paper is organised as follows: In the next section 
we introduce the model of the coupled particle system. Section IIVI is devoted to a Melnikov analysis related to 
homoclinic motion in the system. In the subsequent section we apply a subharmonic Melnikov method to prove the 
persistence of periodic orbits in the coupled dynamics. In section fVTl we present exact results regarding the collective 
particle transport mediated by periodic running solutions. The ensuing current is considered in section fVIIl In IVIIII 
we discuss the synchronisation of the coupled particle system. We prove that for overcritical coupling strength a 
completely synchronised state is attained. Finally we summarise our results and give a brief outlook. 

III. THE COUPLED PARTICLE SYSTEM 

We study the dynamics of two damped, coupled particles evolving in a spatially symmetric (washboard) potential 
of spatial period L = 1 , with equations of motion given by 

qi = -sin(27rgi) - jqi - FsiruTH + 6 ) - n(qi - q 2 ) , (1) 
q 2 = - sin(27rg 2 ) - 792 - Fsin(flt + 9 ) + n(qi - q 2 ) ■ (2) 

An external time-dependent modulation of amplitude F, frequency Q and phase 9o serves for periodic modulations of 
the inclination of the washboard potential. 

We emphasise that the system is unbiased in the sense that the forces averaged over time and space vanish, i.e. 

/dg/J' V^y^U o, i = l,2. (3) 
o o 

where the potential is given by 

V( qi ,q 2 ,t) = — (2 - cos(2tt9i) - cos(2^ 2 )) + F sin(fi t + o )(<Zi + 92) ■ (4) 

Further, notice the exchange symmetry of the system (91,91) <-> (92,92)- 

Before we embark on our study of the coupled particle dynamics we briefly consider the dynamics of a single damped 
and driven particle system, i.e. n = 0. The parameter values, Q — 2.25, F = 1.3, and 9q = 0, are chosen such that 
the dynamics is chaotic. To illustrate the phase flow we utilise a Poincare map using the period of the external force, 
To = 27r/f2, as the stroboscopic time. The system of equations of motion was integrated numerically and omitting a 
transient phase, points were set in the map at times being multiples of the period duration To. The resulting strange 
attractor is shown in the (p = q, q)— plane in Fig. [TJ 

IV. MELNIKOV- ANALYSIS 

In this section we focus our interest on the influence of the coupling between the particles on the character of 
the dynamics. In particular we are interested in a more complex form of unstable motion which may arise in the 
five-dimensional phase space of the system ([U,©, namely hyperchaos. 



q 



Figure 1: Strange attractor for a single particle system, i.e. k = 0. The remaining parameter values are given by SI = 2.25, 
F = 1.3, O = 0, and 7 = 0.1. 



We perform a Melnikov-analysis related to homoclinic motions in system ((J),© for weak coupling, weak forcing, 
and weak damping which is expressed as 



q\ + sin(27nji) = e (-791 - Fsm(flt + 9 ) - n{qi - q 2 )) , 
q 2 +sin(27rg 2 ) = e (~792 - Fsin(f2t + O ) + n(qi - q 2 )) , 



(5) 
(6) 



where the parameter e< 1 indicates the perturbative character of the expressions on the r.h.s. in Eqs. ©,©). For 
e = 0, two unperturbed pendula are obtained 



pi = - sin(27T^) , 
There exist hyperbolic fixed points given by 



Pi , 



1,2. 



{Pi,<li 



0,±; 



1,2. 



Thus viewed in the full p\ — q\ — p 2 — q 2 — space the system has a hyperbolic periodic orbit 

M = (pi,«i,p2,«a,0(t)) = (o,±i,o,±i,nt + o ). 

This hyperbolic periodic orbit is connected to itself by two pairs of homoclinic trajectories given by 

(Pfh (t),qi h (t),P^ h (t),qi h (t),9(t)) = 

7T 



(7) 
(8) 

(9) 



±W-sech(\/27rf) ,±-siri -1 tanh (V^ttj , ±J - secla (V^ttj , 



±— sin 

7T 



tanh ( v 27rt 



, sit + 9 Q 



(10) 



.M has three-dimensional stable and unstable manifolds which coincide along four three-dimensional sets of homoclinic 
orbits, denoted by T^, which can be parametrised as 



{(pth(- T i)^f h (-Ti),pf h (-T 2 ), qf h (-r 2 ), 6> ) 

<E R 1 x T 1 x R 1 x T 1 x T 1 I (ti,t 2 ,0 o ) G R 1 x M 1 x T 1 } 



(11) 



The real parameters ri 2 determine the position on the homoclinic orbits. To determine whether the stable and 
unstable manifolds intersect transversely we compute the Melnikov integrals (for details concerning the Melnikov 



method in higher dimensions see [29]- [32]) 

Mt{rx,r 2 ) 



{Pth(t - n) [-jpf h (t - n) - Fsm(SU + 9 Q ) 



+ « feW (* - n+i) - qf h (t-Ti))]}dt. 



(12) 
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The integration is performed over the homoclinic trajectories given in Eqs. (ITU1) yielding the following integrals 
= J I i^sech (^{t - Ti)) TT^" sech (V^(t - r;)) 



Fsin(Q£ + 6» ) + k f±-sin _1 tanh (y2n(t - 



=F — sin 

7T 



tanh 



(V2tt(* - 7?)) 



tit, i = 1, 2 , r 3 = n, and ^ = 



Using the method of residues we obtain for the case of equal signs 



M^^\ n ,T 2 ) = -2 7A /^ T F sech ( A / - - ) sin(fin) + ^J(At) , 



and 



M^Cn,^) = -2 7 yjT F sech (yf §) Bin(^r 2 ) - ^J(Ar) . 



For unequal signs we obtain 



M 1 (± ^ ) (ti,t 2 ) = -2 7A f, : F *<<•!>( v /|r] sin (fi Tl )__L7(Ar), 



7T fi 



and 



M 2 (± ' T) (r 1; r 2 ) = -2 7A /4 ± F sech ( ) sin (^) + -^J(Ar) 



2 2 



The function I (At) is determined by the integral 



I(At) 



sin -1 [tanh(i + At)] 
cosh(i) 



dt. 



where 



At = V2^(ri - r 2 ) . 
Denoting t = At, the function I(t) has the following properties 



J(0)=0, I(t) = -/(-t) ; 
It can be readily seen that if 



dl(t) 



> , and max | I(t) \= lim | I(t) |= - 

CSC tGR t— >±oo ^ 



7T n 



^> (27\/z3+^ 7 W] cosh U/o o ], 



2 2 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



then there are roots ti,t 2 6 R to the transcendental equations (fT4"l) - (|17l) . Therefore, the Melnikov functions M^ ,ziz ' 
and Af( ±,=F ) have zeros. Furthermore, with 



dM\ 



(±,±) 



9ri 



T-F 1 cos(r2ri) + 7, 



(±.±) 



= —I 



3t 2 

^ = ±F cos(fir 2 )-J, 

OT 2 



(22) 
(23) 



one obtains 



dct DM {± > ±] (ri, r 2 ) = ±[cos(On) - cos(nT 2 )]FI + F 2 cos(firi) cos(Ot 2 ) - P , 



(24) 
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where we used the abbreviations 

F = Osech^--J, and/.--ii. (25) 

Providing that det DM^'^ (r) ^ (det M^'^' (r) ^ 0) then the homoclinic Melnikov vector M {± ' ±s > (M^^) ) has 
simple zeros implying transversal intersections of the stable and unstable manifolds [29]- [32]. 

With regard to the coupling strength, k, the condition in (|21[) imposes also a constraint on the magnitude of the 
product k ■ I(t). Moreover, using the properties in (f2"0"f one obtains the following homoclinic bifurcation threshold 

^=(27\/^ + £]cosh(y|^. (26) 

That is, for F > F c the existence of roots to the transcendental equations (fTl| - ([T7| is guaranteed giving birth to 
homoclinic orbits being of importance for the existence of strange attractors. Note that this bifurcation condition 
limits the maximal value of the coupling strength k but is independent of the value of At. 

With further regard to the influence of the coupling strength k on the behaviour of the stable and unstable manifolds 
the condition for the existence of zeros of the Melnikov functions, as given in (piT]) , confines the magnitude of the 
product k ■ I(At), viz. relates the value of k and At. In fact, taking into account the properties of /(At) (see 
Eq. (|20|) ) one concludes that the larger the coupling strength k the smaller At has to be in order that the value 
of I (At) complies with the condition (1211) . That is, with increasing values of k the points where the stable and 
unstable manifolds intersect each other transversely on the homoclinic orbits get closer to each other, giving evidence 
for the tendency towards synchronisation between the two particle oscillators. Eventually for t\ — Ti the ensuing 
'freezing of the dimensionality' reduces the system to a single driven and damped pendulum (one-and-a-half degree 
of freedom system) which impedes the occurrence of hyperchaos as this is necessarily connected with expansions in 
several directions in phase space simultaneously (several positive Lyapunov exponents). Nevertheless, for sufficiently 
small k, simple zeros of the Melnikov functions with t\ ^ t^ exist allowing for the occurrence of hyperchaos. 

In the limiting case k — ¥ oo in conjunction with At — > the two components of each of the homoclinic Melnikov 
vector degenerate to the Melnikov functions associated with a single pendulum that is damped and periodically driven. 
The corresponding Melnikov function 

M ± (t) = - 7 2W —r i F sech ( J J ^ ) sin(ftT) (27) 



has simple zeros provided the inequality 



2 2 



F 2 y /-|cosh(Vjfi j - CIS) 



2 

is satisfied. Clearly, this type of homoclinic motion cannot represent the source of hyperchaos 



PERIODIC MOTIONS SUBHARMONIC MELNIKOV ANALYSIS 



With view to coherent particle transport mediated by periodic running solutions we investigate the existence of 
periodic orbits associated with collective rotational motion of the particles in the presence of the coupling between 
them where again the perturbative approach is used (see Eqs. , ((6]) ) . For e = each unperturbed subsystem supports 
a one-parameter family of periodic orbits corresponding to rotational motion 



2 1 



TT k 



dn 



/2vr 



t,k 



, i . 

, ± — sin 

7T 



27T 



t,k 



(29) 



where dn and sn are the Jacobi elliptic functions with the elliptic modulus k £ (0, 1) and with period 



(30) 



where K is the complete elliptic integral of the first kind. 
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We investigate periodic motions near the unperturbed resonant tori 7fe 1 .fe 2 , where 

TfcU = {(pi(t + Ti),qt(t + Ti),pi(t + T 2 ),q±(t + T2)),T i e[0,T h ),i = l,2, 
6 = Sit mod (2tt), t £ R} . 

The subharmonic Melnikov functions M-^™ (ri , t 2 ) for Th x ,k 2 are determined by 



(31) 



AC /n (n,r 2 ) 



ro T 



Substituting the expression in (j29l) yields 

/2tt 



Pk ( (* - Ti) [-TP*; (* - - F sin(fii) 

dt, i = 1,2, T3 = Ti , and &3 = fei. 



(32) 



. o ~ /•m T / 



2 F 

7T fej J„ 



m T 



dn 

V2 1 k />m ° T 



TT TT k 



l JO 







^sin" 1 


sn | 



dn 



/27T 



(t - Tj), fc, J dt 
(t — t^, fc^ sin(f2t) dt 
(t - n),ki 



(t — Ti+l), foj+l 



/27T 



(t Ti), hi 



dt. 



where ki satisfies the resonance condition 



/2 27T 
-ki K(fcj) = TOjTo = m,- , i = l,2 
7T ll 



(33) 



(34) 



and mo is the least common multiple of mi and m 2 . Using the Fourier series expressions for the Jacobian elliptic 
functions [28j to evaluate the second integral in Eq. (13U)) one obtains 



Mr /n (r u r 2 ) 



"7 



mo 



V2Tr 3 / 2 k! To] 



niE(fci) 



V2 



FIi(ki,m ,mi,ni = l)sin(firi) + k . I 2 (ti,t 2 ) , 



(35) 



where E(fc) is the complete elliptic integral of the second type. In (|35j) the function Ii(ki, Too, TOi, n\) is non-zero only 
for m — 1 when it is given by 



Ji(fci,mo,mi,ni = 1) = 



27rm ( 



— sech 7rmi 



K(fci) 
■K(fci) 



(36) 



meaning that no periodic oscillations are triggered by the external time-periodic modulation unless n\ = 1. K(fc') 
is the complete elliptic integral of the first kind with k' = \/I — Furthermore, in (|35[) the expression I 2 (ti,t 2 ) is 
determined by the integral 



7 2 (Ar) = dn t, k^j sin" 1 



snl ^l(t + Ar),k 2 



dt. 



(37) 



Denoting i = At, the function I 2 (t), has the following properties 

dl 2 (t) 



h(0) = 0, I 2 (t) = -I 2 (-t) . 



,, >0, h,max{t) = max/ 2 (t) > , and I 2 , m i n (t) = min I 2 (i) < . (38) 



Due to symmetry the second component of the Melnikov vector is determined by 

M m/n (r l ,T 2 ]ki,k 2l m ,m 2l n 2 ) = M[ n/n (r 2 , n ; k 2 , k\ , to , to 2 , n 2 ) 



(39) 
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Furthermore, the Melnikov functions possess the following periodicity 

M™ /n ( n ,7i!) = Mf /n (ri + hT kl , r 2 + l 2 T k2 ) , (40) 

with integers h } fa- 
in particular, if k\ — k 2 = k, mi — m 2 , and n\ = n 2 it holds that t\ — t 2 = lT k leading to a vanishing integral I 2 = 

in Eq. (|37p . Hence, for collective particle motion being frequency locked to the external time-periodic modulation, 

and when, apart from a phase difference, the dynamics of the two particles is identical, there is zero average energy 

exchange between the coupled particles. 

To obtain a criterion for the existence of periodic orbits near the unperturbed resonant tori 7fc 1 ,fc 2 we note that if 

r, ( 4 tylq , . V2 T , A fifciK(Jfci) , / K(jfei)\ /t . 

F> {■<7^^ niEik > ) + K ^ h ™' m ) ^7 cosech (' m, skl) ' (41) 

then there exist roots n, t% £ [0, 27r/f2) to the transcendental equations ([35]) and (|39|). Hence, M m l n — (AI™^ n , M™^ n ) 
has zeros (ti,t 2 ). Furthermore, as DM m /"(ri, t 2 ) ^ 0, the system CO),© possesses subharmonic orbits of order mo 
near each of the unperturbed resonant tori 7fc 1 .fc 2 [29]- [32]. 



VI. FEATURES OF TRANSPORT 



In the previous section we proved for weak coupling, weak forcing, and weak damping the existence of periodic 
motion forced by the external time-periodic modulation facilitating the subharmonic Melnikov method. In particular, 
it turns out that periodic orbits are characterised by T — j To, with integer j, and T is the period of rotational 
motion of the particles and T = 2~k /VL is the period of the external time-dependent modulation. Departing from 
such a perturbational approach we now investigate the features of running solutions accomplishing coherent particle 
transport. Periodic running solutions are characterized by 

q i (t + T)=q i (t) + m i , q^t + T) = q t (t) , i = 1,2, (42) 

with some period T and non-vanishing average velocity 

(4i) = ^J dtqi(t) = Ci^0, (43) 

where a constant Ci is positive (negative) when particle i runs to the right (left), that is m., > (rrij < 0) and i = 1,2. 

As noted in section HV1 for in-phase motion qi(t) — q 2 (t), that is for identical dynamics of the two particles, the 
coupling between the particles vanishes and the dynamics is described by the system of a single particle. Crucially, for 
our choice of the values of the parameters the single-particle system possesses a strange attractor. Therefore, in-phase 
running solutions supported by stable attractors associated with periodic transporting motion are excluded. 

As the character of other possible running solutions is concerned we state the following 

Theorem: For the coupled system given in Eqs. ((J),© with k > and 7 > 0, for the periodic running solutions 
determined by Eqs. (g^J) and (g3J it holds that 

(i) During one period, T, the two particles run over an equal distance, that is raj = m 2 . 

(ii) No (nontrivial) periodic motion is possible in the absence of the time-periodic external modulations, i.e. if F(t) = 0. 

(iii) For solutions, being frequency-locked to the external time-periodic modulation, F(t) = F sm(flt + 6o), with period 
T = 2tt/0, the distance between the particles performs periodic oscillations due to 

qi(t + T)-q 2 (t + T)= qi (t) - q 2 (t) , (44) 

and the period T is determined by 

T = 2lT 0l (45) 

with integer I > 1. 

(iv) The coordinates obey 

9i (t + \t\ = q 2 (t) + k, q2 (t+h\= qi (t) + k , (46) 



with some integer k ^ 0, and, hence 



q. l (t + T) = q l {t)+2k, i = l,2. 



(47) 



Proof: Suppose that there exists a period T such that the system ([lj.© has a running solution of the form given in 
(|42l) . Multiplying Eq. ([1]) by q± and Eq. <j2j> by q 2 and adding the two resulting equations we obtain 



d 



111 1 K 

2«i + - ^ cos ( 27r '? 1 ) _ 2^ cos ( 27r 9i) + 2 ( 9l ~ 92 ) 



= -F(t)(<h+<b)-7(«l+«£), 
where F(t) = F(t + To) = Tsin(fii + 9q) and To = 2-k/Q,. Integrating over one period, T, yields 



(48) 



nT 

i 

(n-l)T 
nT 



111 1 K 

2^ + ~$ - 27 cos(27r gi ) - — cos(27rgi) + -{q x - q 2 f 



dt[F{t){qi+q 2 )+7{ql + ql)] . 

/(n-l)T 

With the running solution as given in (|4"2l the expression in (|4"9")l reduces to 



'(n-l)T 

We denote the r.h.s. of Eq. ([50)) by 



2(9i - 92) 2 



(n-l)T 



dt [F(t)(?!(t) + q 2 (t)) + 7 (#(t) + q 2 2 (t))] 



nT 



In = - I dt [F(t)( qi (t) + q 2 (t)) + 7 (<Z 2 W + q 2 2 (t))] 

/(n-l)T 



(49) 



(50) 



(51) 



Taking into account the property in Eq. (|4"5|) , sinusoidal F(t), and discarding the second negative-definite term 
under the integral it follows that there exists a constant C > such that the value for /„ can be bounded from above 
as follows: 



In < FT(\ d I + I C 2 I) = C = const. 
On the other hand, for the l.h.s. of Eq. (|50[) we derive 

iT 



(52) 



(n-l)T 



Hqi - (12)' 



(91 (nT) - 92 (nT)) 2 - {q x {{n - 1)T) - g 2 ((n - 1)T)) 2 



= ^{[(9i(0)+nmi)- fe(0)+nm 2 )] 2 

- [((71 (0) + (n- l)m x ) - (<h(0) + (n - l)m 2 )] 2 } . 

Using the notation Ag(0) = qi(O) — 92(0) and Am = mi — to 2 one arrives eventually at 

A K)Ti = [2Ag(0) + (2n - l)Am]Am . 



(53) 



(54) 



For Am ^ the term A Kjn grows monotonically with increasing n and therefore, A K .„, cannot be balanced by the 
corresponding finite counterpart /„ appearing on the r.h.s. in Eq. (|50[) . A Kj „ attains a finite value, i.e. becomes 
?i-independent only, if Am = 0, i.e. mi = m 2 which concludes the proof of statement (i) in the above Theorem. 



For the proof of part (ii) we note that Am = in Eq. (|5"4"|) implies 

A re ,„ = / d -(qi - q 2 y 

J(n-1)T 11 



= 0. 



(55) 
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so that due to Eq. (|50|) the following relation is true 



dt [F(t)( qi (t) + & (*)) + 7(9? (*) + 9 2 2 (*))] = o , 



-1)T 

which in the absence of the external time-dependent modulation, i.e. F(t) = 0, reduces to the condition 



dt(qf(t) + q^t)) =0. 



(56) 



(57) 



(n-l)T 



This condition can only be satisfied if q\ = q 2 = 0, q\ = const., and q2 = const, which, however, contradicts the 
assumption of (nontrivial) periodic solutions. Thus, to obtain periodic solutions F{t) cannot be zero which proves 
part (ii) of the Theorem. 



In order to prove the remaining statements (iii) and (iv) we introduce the following variables 

x = qi-q 2 , y = qi + q2, 



(58) 



for the difference and sum of the coordinates q\ and q 2 respectively. The corresponding equations of motion expressed 
in the new variables read as 



x = — 2kx — jx — 2 sin(7rx) cos(7ry) = f(x, x, y) 
V = -jy-'2F(t)-2cos(irx)sm(iry) = g(x,y,y) 



(59) 
(60) 



While only the second equation (I60|) contains the time-periodic external modulation, the impact of the coupling, 
related with the parameter k, enters only the first equation (|59[) where the corresponding term, — 2kx. serves as a 
harmonic restoring force that keeps the motion of the difference variable x bounded. In contrast, the sum variable 
y can undergo ongoing rotational motion triggered by the external modulation of period T , viz. F(t) = F(t + Tq). 
Accordingly, periodic solutions y(t) of Eq. (|6"0|) are supposed to be frequency-locked to (multiples of) the external 
periodic modulation F(t), and using (|42[) and the proven result from (i), i.e. mi = ni2 = m, the periodic running 
solutions of the difference variable y(t) attain the form 



y(t + IT )= y{t) + 2m , y(t + IT ) = y(t) , and I > 1 



(61) 



Given further the reflection symmetry of the difference equation (|59p . x <-> — x, symmetric (limit cycle) oscillations 
of x(t) around x = of some period T with 



x(t + T) = x(t) , x(t + T) = x(t) . 



are possible for which the following holds 



T 



x \ t + ^) =-*(*). * *+T =-i(*) 



T 



(62) 



(63) 



In order that Eq. (|59|) supports solutions obeying the relation (|63|) its r.h.s., f(x, x, y), needs to satisfy the reflection 
symmetry 



f x t + 



T 



,x[t + 



T 



y\t + 



T 



= -f(x(t),x(t),y(t)). 



Taking into account the r.h.s. of Eq. (|59|) this results in the following condition 

T 



cos 



ny[t + 



= cos [iry (*))] 



From Eq. ([BT]) we obtain y(t) = y(t + I Tq) — 2m which, after substitution in the r.h.s. of Eq. (|65|) gives 

T 



cos 



ny [t + 



= cos [iry(t + IT )} 



(64) 



(65) 



(66) 



from which follows the relation between the frequencies of the periodic motions in x(t) and y(t): 

T = 21 T . 



(67) 
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Using the relation cos[7ry(t + 21 To)] = cos[7r(y(t) + 2 • 2m)] = cos[7ry(i)] one readily verifies the validity of the 
condition 

/ (x(t + T),x(t + T), y(t + T)) = f(x(t), x(t),y(t)) , (68) 

that has to be satisfied further by the r.h.s of Eq. (JSHJ) in order to support periodic solutions which comply with 
Eq. (EU). 



Furthermore, given the relations in (f6"2")l , (|6"5|) and (f6"T)) it holds that 

x(t + lT ) = (-l) l x(t), x(t + lT ) = (-l) l x(t) (69) 

for integer I. It remains to verify that the expressions (|5T|) and (|6"9")) for periodic solutions x(t), y(t), that are frequency- 
locked to the external time-periodic modulation F(t), leave the r.h.s. of the sum equation (|60[) invariant. We obtain 

g(x(t + IT ), y(t + I T ),y{t + IT )) = -jy(t + IT )- F(t + lT )-2 cos[7ra(t + I T )] sin[^(< + I T )] 

= -iy(t) - F(t) - 2cos[7r(-l) J x(t)] aia[w(y(t) + 2m)] 
= -iy(t) - F(t) - 2cos[7ra;(t)] sm[iry(t)] 

= g(x(t),y(t),y(t)), (70) 
ensuring invariance which completes the proof of the statement in (iii) . 

Finally, combining the expressions in (|5T|) and (|6"9")) one infers that in original coordinates the following holds 

qi (t + lT ) = q 2 (t)+m (71) 
q 2 (t + lT ) = <&(*)+ m (72) 
q t {t + 2lT {) ) = qi (t) + 2m, z = l,2, (73) 

which concludes the proof of part (iv) of the Theorem. □ 

Corollary: No counterpropagating running states exist as mi and m,2 can only have the same sign implying that the 
two particles run in the same direction. 

In other words, the theorem states that the periodic transport scenario is characterised by coordinated energy 
exchange between the two particles proceeding such that they alternatingly surmount the barriers of the periodic 
potential. We illustrate the features of coherent transport in the periodic regime in Fig. [2] displaying the time 
evolution of the coordinates 31,2^) of the two particles for coupling strength k — 0.46, that is within the window of 
periodic motion of the bifurcation diagram (see Fig. 2]). To illustrate the frequency- locking between the coordinates 
and the the time-periodic external modulation, F(t) — Fsin(Clt), the latter, oscillating around q = 1564 with unit 
amplitude F = 1, is also shown in Fig.O Fig. [3] displays the evolution of the corresponding trajectory in the spatially 
periodic potential landscape U(q\,q-2) = [2 — cos(27r<7i) — cos(27r<7 2 )]/(27r). Moreover, the temporal behaviour of the 
coordinates follows the relations given in Eqs. (|44I) - (|4"T|) . 

In the following we deal with the existence, uniqueness and approximation of periodic solutions given in (|42|) where 
m i — m 2 = m - To this end we cast the system ([I},© in matrix notation 

Z(t)=AZ(t)+B{t,Z(t)), (74) 

where Z = (pi,P2j <Zi, <?2) T and the matrix A is given by 



A = 



V 



7 





— K 


K 





-7 


K 


— K 


1 














1 









The inhomogeneity is given by 

B(t, Z{t)) = (- sin(2^i(t)) - F(t), - sm(2nq 2 (t)) - F(t), 0, 0) T . (75) 
The solution to ([74")) with initial conditions Z(0) reads as 

Z(t) = exp(At)Z(Q) + [ dr exp[A(t - t)]B(t, Z{t)) . (76) 
Jo 
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5020 



Figure 2: Time evolution of the coordinates qi,2(t) of the two particles representing a periodic running state being frequency- 
locked to the external time-periodic modulation. The parameter values are given by = 2.25, F = 1.3, 9o — 0, 7 = 0.1, and 
K — 0.46. For comparison, sin(ttt) oscillating around q — 1564 with unit amplitude F = 1 and frequency Q = 2.25 is shown 
(dashed line). 



1567.6 




1559.6 



1567.2 



1559.9 



Figure 3: Evolution of the trajectory associated with the two particles in a periodic running states in the spatially periodic 
potential landscape U (91,92)- The parameter values are the same as in Fig. [2] 
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Periodic solutions of period T, as given in (|42p with mi = — m, fulfil the relation 

Z(T) = Z(0) + M, (77) 
where M = (0,0, m 1 m) T . Substituting Z(0) = Z(T) — M into (|76f yields a nonlinear integral equation for Z(T): 

i-T 



Z(T) = (I-exp(AT))- 1 



■exp(AT)M + / dreoq)[A(T-r)]B(T,Z(r)) 
Jo 



(78) 



Concerning the solutions Z(t) of the nonlinear integral equation (l78l) we state the following 

Theorem: The integral equation (|78p with T 6 [0, f] possesses a unique solution. Moreover, if a sequence of functions 
{Z„} is defined inductively by choosing any Z £ C([0,t];R 4 ) and setting 



Z n+ i = (I-exp(AT)y 1 



-exp(AT)M + / drexp[A(T-r)]B(T,Z n (r)) 
Jo 



(79) 



the sequence {Z n } converges uniformly on / := [0,t] to the unique solution Z of (J75J). 

Proof: The proof utilises Banach's fixed point theorem. We consider the space of continuous functions Z from / to 
R 4 equipped with an exponentially weighted metric (the Bielecki metric [33| ) 



- II Z(T)-Z(T) || 

d a (Z,Z):= sup — — 

Te[o,t] exp(aT) 



(80) 



where || • || denotes the Euclidean norm on R 4 and a > is a constant. The corresponding norm is determined by 



Z \\ a = SUp 



Z(T) 



Te[o,t] exp(aT) ' 

C([0, i]; M 4 , || • || Q ) is a Banach space. A solution to (|78|) is a continuous function Z : I H> R 4 . Let S* 



5 : C*([0,t];R 4 ) h> C([0,t];R 4 ) 



be defined by 



[5Z](T) = (I-exp(AT))- 1 



-exp(AT)Af+ / drexp[A(T- r)]B(r,Z(r)) 
Jo 



(81) 



(82) 



(83) 



Then the fixed point of S represents the solution to (|78[) . Hence we must prove that there exist a unique Z such that 
SZ = Z. 
Let 



= sup || exp(A(T - r)) || , K 2 = sup || (I - exp(A(T)))- 1 || , 

(T,r)e[0,t]x[0,t] Te[0,t] 

with a suitable matrix norm || • ||. Further, using 

|sin(27ra;) — sin(27rj/)| = 2|cos[7r(x + y)} sin[7r(a; — y)]\ 
< 2|sin[7r(x-y)]| < 2n\x-y\, 

implies that 

\\B(T,Z(T))-B(T,Z(T))\\<2tt\\Z(T)-Z(T)\\, VTe[0,t], (Z,Z)eR 4 , 
and B(T, Z{T)) is defined in (|75|) . Let a = 2ttKiK 2 S, where 5 > 1 is an arbitrary constant. 



(84) 



(85) 
(86) 
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Figure 4: (Colour online) Bifurcation diagram as a function of the coupling strength and remaining parameter values: Q = 2.25, 
F = 1.3, O = 0, and 7 = 0.1. 



In order to apply Banach's fixed point theorem we show that S is a contractive map with contraction constant 
1/6 < 1. For any Z, Z € C([0,i];R 4 ) one gets 

^sz,si) - sup ll[^]m-[^(r)ll 



re[o,t] exp(aT) 



< K X K 2 sup ]— f T dr\\B(T,Z(r))-B(T,Z(T))\\ 

TeF0.il exp(aT) / 



Te[o,i] exp(aT) J 

1 /■ T ^_^J^)-^)H 



< 27rJfi^ 2 sup ^— / dr || Z(r) ~ Z{t) 

Te[o,i] exp(aT) 7 



1 r 

2itK\K2 sup / drexp(aT) 

TeF0.il exp(aT) 7 



Te[o,t] exp(aT) J 1 ' exp(aT) 

1 f T 

< 2nK 1 K 2 d a (Z, Z) sup — — / drexp(ar) 

Te[o,t] exp(aT) J 

= a Q (^, Z) sup [1 — exp(— al )\ = - 

a Te[o,t] V a d 

< I [l-cxp(- at)} d a (Z,Z)<jd a (Z,Z). (87) 



Since S > 1 one concludes that 5 is a contractive map so that Banach's fixed point theorem applies assuring the exis- 
tence of a unique fixed point Z of S. Moreover, from Banach's theorem follows that the sequence {Z n }, defined in (|7^|) 
converges uniformly in the norm || • || a to that fixed point Z. □ 



VII. EMERGENCE OF A CURRENT 



The dynamics of the coupled particles exhibits very rich and complex behaviour and depending on the value of 
the coupling strength k > 0, one finds periodic, or aperiodic (quasiperiodic and/or chaotic) solutions in the long 
time limit. The character of the phase flow evolving in a five-dimensional phase space is conveniently displayed by a 
Poincare map using the period of the external force, To = 27r/f2, as the stroboscopic time. The system of equations of 
motion was integrated numerically and omitting a transient phase points were set in the map at times being multiples 
of the period duration To. In Fig. 3] the bifurcation diagram as a function of the coupling strength is depicted. 

Particle transport is quantitatively assessed by the mean velocity, v m , which we define as the time average of the 
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Figure 5: (Colour online) Mean velocity v m /vo as a function of the coupling strength and the remaining parameter values are 
as in Fig. [4] 

ensemble averaged velocity, i.e. 

v m = ^- ( dt' <«i,„(t') + «a, n (0> » ( 88 ) 

Is Jo 

with simulation time T s and with the ensemble average given by 

1 N 

(Mt)) = -Y,vn(t). (89) 

n=l 

Here N denotes the number of particles constituting the ensemble with associated random initial conditions q n (0) 
and p n (0) that are uniformly distributed over the period of the potential. We express v m in terms of the ratio of the 
spatial and temporal periods L/Tq = vq with vq — 0.358 being the velocity for running solutions that advance by one 
spatial period during one period duration of the external field. 

We recall that in the uncoupled case, n = 0, the dynamics is characterised by a strange attractor. Concerning 
transport, the associated scenario is that of chaotic running solutions which however, averaged over time and ensemble, 
do not contribute to a non-zero current. 

For increasing coupling strength, k > 0, two typical scenarios arise: pinned and running states. In the former state 
the motion proceeds at most over a finite number of spatial periods whereas in the latter state motion is directed 
and unrestricted in the spatial dimension. In terms of the phase flow, running asymptotic solutions correspond to 
phase locked periodic attractors transporting a particle with velocity v — m/n over m spatial periods of the potential 
during n period durations To of the external periodic field. Running asymptotic solutions may also be supported by 
aperiodic attractors. 

In the bifurcation diagram associated with the dynamics shown in Fig. U one recognises vertically extended stripes 
covered densely with points corresponding to non-phase locked aperiodic attractors and several periodic windows as 
well as period-doubling cascades to chaos. These features of the phase flow are readily attributed to the resulting 
mean velocity of the net motion (depicted in Fig. [5]). The ensemble average is taken over an ensemble of N = 5000 
trajectories with uniformly distributed initial conditions q(0) and q(0) = v(0). For computation of the long-time 
average the simulation time interval for each trajectory is taken as T s = 5 x 10 5 ~ 1.8 x 10 5 x To. We notice almost 
in the entire K-range vanishing mean velocity v m — 0. The exception is the interval 0.418 < k < 0.504 for which the 
solutions are associated with multiple coexisting attractors lying in a fairly extended periodic window in Fig. [4] (In 
the periodic window in the interval 0.082 < n < 0.094 all solutions are pinned on non-transporting attractors, and thus 
the resulting current is zero.) We observe that the current oscillates wildly with changes of its sign. That is, oppositely 
running solutions attributed to branches of the various coexisting periodic attractors contribute to the mean velocity 
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with different weight with the effect that either those with positive, v = vq, or negative velocity, v — —vq, dominate 
yielding the window of changing mean velocity v m . For increasing re > 0.504 these periodic attractors are destroyed 
by way of crisis after passage through a period-doubling route to chaos. 

The directed motion results from a lowering of the dynamical symmetry caused by the external modulation field 
[16], [18]. That is, even though the potential and the external modulation field are with respect to space and time 
symmetric respectively, with the choice of a fixed phase 9o the symmetry of the flow is reduced and a phase-dependent 
net motion is found. (Note that additional averaging over the phase 9q yields zero mean velocity.) Due to symmetry 
reasons it holds that the sign of the mean velocity is reversed upon the changes 9q = — > 9q — t and F — > — F 
respectively. However, there exists a phase < 9q < n for which symmetry between the two coexisting periodic 
attractors supporting solutions with velocities of opposite sign, v and — vq, is restored and therefore the net motion 
vanishes. 

VIII. COMPLETE SYNCHRONISATION 

In this section we consider the synchronisation features of the coupled system which are related with the properties 
of the difference system given in Eq. (]59| . 

Using the notation X — (x,±) T and F = (0, — 2 cos[7ry] sin[7ra:]) T the solution to Eq. (|59p with initial conditions 
X(0) is expressed in form of an integral equation 

X(t) = exp(At) Jf (0) + [ dr exp[A(i - t)]F(t, X(t)), (90) 
Jo 

where cos[7ry(i)], contained in the inhomogeneity F, acts as a time-dependent driving term. The matrix A reads as 

1 

~2k -7 



The eigenvalues of A are given by 



A 



± 



I (iT Vj 2 -Sk) , (92) 



and fulfil the inequality Re(A±) < —T, T > with T = (7 + ^Ji 1 - 8k)/2 if k < 7/8 and T = 7/2 if re > 7/8. Hence, 
there exist a K > such that the relation 

|| exp(Ai) \\<K- exp(-H) (t > 0) (93) 

holds. || • || is a suitable matrix norm. Further, note that due to 

|cos[7ry] sin[7rx] | < |sin[7rx]| < 7r|x| , (94) 

the inhomogeneity F satisfies 

\\F(t,X(t))\\<2n\\X(t)l (95) 

with || ■ || denoting the Euclidean norm on R 2 . 
First, we prove the following theorem. 

Theorem: Consider the integral equation (|9"0"]) . If 

T>2ttK, (96) 

then for all initial conditions X(0) the solutions X(t) to ([9H| approach zero as t — >• 00. 
Proof: Concerning the norm of X one gets the following bound 

|| X(t) || < Xexp(-rt) || A(0) || +K f drexp(-r(t-r)) || F(t,X(t)) \\ 

Jo 

< Kexp(-Tt) || X(0) || +2wK [ drexp(-r(t - r)) || X(r) || (97) 

Jo 
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Setting G(t) = exp(rt) \\X(t)\\ leads to 

G(t) < K || X(0) || +2ttK [ (1tG(t). (98) 

Jo 



Using Gronwall's inequality we obtain 



2irK / dr 



= K\\X(0)\\exp[2nKi\. (99) 



G(t) < K || X(0) || exp 
Therefore, 

\\X(t)\\<K\\X(Q)\\exp[-(T-2irK)t] , (100) 

and as r > 2irK it follows that 

hm || X(t) ||= 0. (101) 

□ 

Corollary: As -X'(i) = (<c(i), x(t)) = (qi(t) — 92 (*) ! 9i (*) — 92 (*)) one concludes that under the condition (IM|) for 
the system of equations (HJ),© the motion completely synchronises asymptotically and exponentially fast regardless 
of the initial conditions. 

We remark that nothing can be said about the character of the completely synchronised state of the system which 
can be regular but also chaotic. 

A. Transverse stability of the synchronous attractor 

The exchange symmetry properties of the system ((J),© implies that the hyperspace q\ = 92 and q\ = 92 constitutes 
the synchronisation manifold. The latter is invariant, viz. if at any moment of time the coordinates and velocities of 
the two particles coincide, then they will do so for all times. Within the synchronisation manifold the dynamics is 
governed by the equation for a single particle. 

As established in section Mil for our choice of parameter values the dynamics of a single oscillator is characterised 
by a strange attractor. 

In the following we investigate the stability of the synchronisation manifold. To this end we utilise the difference 
system given in Eq. (|59j) in terms of which the synchronisation manifold is represented by the origin (x, x) = (0, 0) 
in phase space. Therefore, stabilisation of the origin in the difference system is associated with complete (globally 
identical) synchronisation of the motion of the two particles. 

The stability of perturbations transverse to the synchronisation manifold is governed by the largest (transverse) 
Lyapunov exponent associated with the difference system [34| . To be precise we study the corresponding tangent 
equation 

Sx = —2k5x — 7<5± — 2tt cos(TTc(t))dx , (102) 

where c(t) = qi(t) = 92 (t) corresponds to completely synchronous motion. In Fig.[f)]the largest (transverse) Lyapunov 
exponent as a function of the coupling strength is shown. The largest transverse Lyapunov exponent becomes negative 
for K > 5 indicating complete chaotic synchronisation for sufficiently strong coupling strength in compliance with the 
effect of 'freezing of the dimensionality' occurring for increasing value of k (see Section HV|) . Notably, as in the 
completely synchronised system the coupling between the two particles vanishes, the occurrence of hyperchaos is 
impossible. Furthermore, although the dynamics in the coherent transport regime, appearing for coupling strengths 
in the range 0.418 < k < 0.504, is regular, the corresponding transverse Lyapunov exponents are positive as the 
dynamics does not evolve in the synchronisation manifold and the distance between the particle rather changes 
periodically according to Eq. (|4"4")l . 



IX. SUMMARY 



We have considered the nonlinear dynamics of two coupled, damped and time-periodically driven particles evolving 
in a symmetric and spatially periodic potential. One aspect of our study has concerned the emergence of unstable, 
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Figure 6: The largest transverse Lyapunov exponent in dependence of the coupling strength k. Stability of the synchronised 
chaotic state is assured for k > 5 for which A t becomes negative. The remaining parameter values are: = 2.25, F = 1.3, 
0o = 0, and 7 = 0.1. 

irregular motion in the coupled dynamics in the five-dimensional phase space. In order to show the existence of 
transverse homoclinic orbits we have invoked a homoclinic Melnikov analysis. The latter has enabled us also to predict 
the occurrence of homoclinic bifurcations giving birth to homoclinic orbits being of importance for the existence 
of strange attractors. Further, we have paid special attention to the influence of the coupling strength on the 
synchronisation features of the coupled particle system. For relatively weak coupling strengths the dynamics in the 
five-dimensional phase space may exhibit simultaneous expansions in several directions leading to hyperchaos. On 
the other hand, we have proved that for overcritically strong coupling strengths the motion of the two particles 
completely synchronises. The resulting 'freezing of the dimensionality' reduces the system to a single driven and 
damped pendulum, i.e. a one-and-a-half degree of freedom, impeding the emergence of hyperchaos. 

Another aspect of our study has dealt with directed particle transport provided by running periodic solutions. 
Notably, for parameter values for which each uncoupled particle exhibits chaotic dynamics, nevertheless, varying the 
coupling strength, a fairly extended window of regular periodic motion appears in the bifurcation diagram. As the 
periodic motion within this window of regular motion is associated with multiple coexisting transporting attractors, 
contributing with different weight to the net flow a non- vanishing current arises (as a result of associated symmetry- 
lowering due to the choice of the phase of the external modulation). We have given exact results concerning the 
features of collective periodic transporting solutions being frequency-locked to the external time-periodic driving. In 
detail, we have shown that coordinated energy exchange between the particles proceeds such that they are enabled to 
overcome - one particle followed by the other - consecutive barriers of the washboard potential resulting in collective 
directed motion. 

Finally, we remark that it is certainly of interest to extend the present study beyond the case of two coupled 
particles to infer general features of coherent many-particle transport in periodic landscapes. 
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